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1.  Introduction 


Ice  is  the  most  well-studied  mineral/™5  The  water  molecule  by  itself  is  capable  of  many  electronic, 
vibrational,  and  rotational  absorption  features.  In  liquid  and  solid  form,  the  rotational  transitions  are 
modified  by  interactions  with  nearby  molecules,  and  the  vibrational  bands  are  shifted  in  wavelength. 
In  solid  form,  ice  shows  additional  structure  due  to  lattice  effects  in  all  crystalline  and  amorphous 
states.  All  transitions  are  functions  of  temperature. 

Figure  1  shows  the  absorption  coefficient  (defined  in  Section  2)  of  ice  and  water  near  water’s  freezing 
temperature  (273K)  from  the  ultraviolet  to  the  far  infrared.6  Both  are  generally  similar  because  they 
are  dominated  by  electronic  and  vibrational  absorptions  of  the  water  molecule  itself.  Only  longward 
of  about  8  pm  do  the  two  curves  differ  significantly.  These  differences  are  the  result  of  solid-state 
lattice  vibrations  in  ice  that  do  not  occur  in  water.  Water  and  ice  are  most  transparent  in  the  visible 
portion  of  the  spectrum,  where  the  absorption  coefficients  reach  a  minimum  around  0.45  pm.  This  is 
the  region  where  so  much  attention  has  been  focused  on  “blue-green”  lasers  for  underwater 
communication. 


Figure  1 .  Absorption  coefficient  a(K)  of  water  and  ice  from  the  ultraviolet  to  the  far  IR.  Ice  data  are 
from  Warren,6  and  water  data  are  from  Hale  and  Querry.7 
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At  temperatures  below  about  137K  in  a  vacuum,  ice’s  lattice  structure  is  amorphous.  Between  137 
and  197K,  ice  is  cubic,  and  above  197K  it  is  hexagonal.  Any  feature  that  originates  in  the  water 
molecule  (like  the  OH  stretch  at  3.16  pm)  will  also  occur  in  the  spectrum  of  ice,  though  there  will  be 
minor  differences  as  a  result  of  lattice  interactions.  Any  feature  that  originates  within  the  ice  lattice 
itself  (like  the  13-  and  45-pm  features)  will  be  absent  from  water’s  spectrum. 

In  general,  the  colder  the  ice,  the  lower  the  absorption  coefficient.  Figure  2  shows  the  absorption 
coefficients  for  warm6  and  cold  ice8  (250  and  80K,  respectively).  Although  the  overall  structures  are 
the  same,  there  are  significant  differences  in  the  spectra,  and,  in  all  cases,  cold  ice  is  less  absorbing 
than  warm  ice. 

In  the  2-20  pm  range,  the  spectra  are  dominated  by  the  four  features,  and  no  others  are  possible  for 
pure  ice.  The  large,  narrow  peak  at  3.16  pm  is  due  to  O-H-O  symmetric  stretch  vibration  (vj).  The 
6.1 -pm  feature  is  due  to  the  V2  bending  mode,  and  the  feature  at  4.6  pm  is  a  combination  of  V2  and  a 
weak  libration  mode.  The  broad  peak  centered  at  13  pm  is  a  lattice  absorption.  There  are  additional 
weak  absorption  features  shortward  of  2.5  pm  (Figure  1),  but  they  normally  are  not  important  in  the 
AI&T  environment. 


Figure  2.  Absorption  coefficient  cx(k)  of  warm  (250K)6  and  cold  (80K)8  ice  based  on  tabulations  of 
Warren6  and  Hudgins  et  al.8 
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2.  Transmission,  Absorption  Coefficient,  and  Thickness 


To  compute  transmission,  t,  the  absorption  coefficient,  a(L),  and  the  thickness,  D,  must  be  known: 

t  =  exp(ADa)  =  exp(-x)  (1) 

where  the  product  Da  is  called  the  optical  depth  or  optical  thickness  x.  D  has  units  of  distance,  and  a 
has  units  of  1 /distance,  typically  cm  and  cm"1,  respectively.  Clearly,  x  is  unitless.  Here,  as  in  all 
subsequent  calculations,  we  will  assume  normal  incidence  and  ignore  the  effects  of  reflections  from 
the  free  surfaces  of  the  ice  layer  (a  few  percent  according  to  Fresnel’s  equations).  The  subject  of 
surface  reflections  is  discussed  in  Section  8. 

The  transmitted  irradiance  I  (colloquially  “brightness,”  “intensity,”  “power,”  etc.)  is  given  by 

I  =  Io  t  =  Io  exp(-Da)  =  Io  exp(-x)  (2) 

where  Io  is  the  incident  irradiance.  According  to  Lambert’s  Law,  the  absorption  coefficient  is  related 
to  the  imaginary  part  of  the  complex  index  of  refraction,  N  =  n  +  ik,  by 

a  =  47rkA,  (3) 


The  presence  of  X  in  Eq.  (3)  means  that  a (X)  decreases  with  increasing  wavelength,  other  things 
being  equal.  Note  that  water’s  absorption  coefficient  (Figure  1)  drops  off  roughly  as  Ifk  longward  of 
about  10  pm  as  Eq.  (3)  suggests  because  there  are  no  discrete  absorptions  in  water  at  longer 
wavelengths.  Note  that  Eqs.  (1),  (2),  and  (3)  are  generic, and  apply  to  virtually  any  solid  or  liquid. 

Equation  (1)  can  be  used  in  a  number  of  ways.  If  a (A,)  is  known — say  from  the  optical 
constants — then  D  can  be  determined  from  the  transmission.  Similarly,  the  transmission  can  be 
computed  for  any  given  thickness.  If  D  and  t  are  known,  a(k)  can  be  retrieved  and  used  for  future 
predictions. 

The  quantity  1/a  is  called  the  skin  depth,  (D,  and  corresponds  to  an  optical  depth,  x,  of  unity.  To 
quickly  estimate  how  thick  a  layer  of  ice  is  necessary  to  produce  significant  absorption,  one  need  only 
consider  the  skin  depth  (Figure  3).  For  example,  in  the  blue-green  region,  an  optical  depth  of  1 
occurs  in  a  slab  that  is  about  3300  cm  thick.  At  the  3.1-pm  absorption  peak,  however,  only  about  0.3 
pm  of  material  is  needed. 
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3.  Determining  the  Absorption  Coefficient,  a(X) 


The  absorption  coefficient  is  usually  computed  from  the  known  optical  constants  (Eq.  3),  or  the 
measured  transmission  (Eq.  1). 

The  term  “optical  constants”  refers  to  the  spectrum  of  complex  indices  of  refraction  N(E)  =  n(>.)  + 
ik(X-),  though  it  is  sometimes  applied  to  (and  compiled  as)  the  complex  dielectric  function  e  =  e’  +  ie”. 
The  two  quantities  are  related: 


(n  +  ik)2  =  e’  +  ie”, 


(4) 


which  is  easily  solved  for 


e’  =  n2-k2  and  e”  =  2nk 


or  equivalently 

n  =  (1/V2)  [(e’2  +  e”2)1Q  +  e’]1/2  &  k  =  (1/V2)  [(e’2  +  e”2)1/2-  e’]m 
Broadly  speaking,  n  controls  reflection,  refraction,  and  scattering,  and  k  controls  absorption. 

In  the  discussion  that  follows,  we  will  concentrate  on  ice’s  spectrum  in  the  1-20  pm  region  where 
many  space  sensors  operate. 


(5) 


(6) 
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4.  Cryogenic  Ice 


Ice  formed  at  extremely  low  temperatures  is  usually  amorphous  and  probably  does  not  occur  naturally 
on  Earth.  The  coldest  temperatures  on  Earth  occur  at  the  summer  polar  mesopause  (-150K),  where 
rare  ice  clouds  called  noctilucent  clouds  form.9  The  clouds  are  thought  to  be  composed  of  hexagonal 
ice  or  mixtures  of  water  and  sulfuric  acid.  Amorphous  ice  may  have  formed  naturally  in  the  outer 
solar  system  or  in  the  cores  of  dense  molecular  clouds  in  space.  Naturally  occurring  ice  clouds  such 
as  cirrus  are  composed  of  hexagonal  ice.10 

Cryogenic  ice  is  amorphous,8,11  and  its  absorption  coefficient  is  temperature  dependent.8  Figure  4 
shows  the  transmission  spectrum  of  a  0.13-pm-thick  film  of  ice  deposited  at  1  OK  in  a  vacuum. 


Wavelength  (/*m) 

Figure  4.  Transmission  spectrum  of  a  0.13-;m-thick  film  of  ice  deposited  in  vacuum  onto  a  10K 
surface  based  on  tabulations  of  the  Jena  group  (vvww.astro.uni-jena.de). 
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The  temperature  dependence  of  amorphous  ice’s  absorption  characteristics  have  been  well  studied. 
The  results  generally  agree;  though,  perhaps  owing  to  different  experimental  conditions  and  different 
types  of  amorphous  structure,  there  are  some  variations.  Shortward  of  about  4  pm,  the  absorption 
coefficient  increases  with  increasing  temperature,  typically  by  about  25%  between  10K  and  140K. 

Figure  5  shows  transmission  spectra  of  80K  amorphous  ice  for  a  variety  of  thicknesses  in  the  2.5-12 
pm  region  based  on  the  work  of  Hudgins  et  al.8 

Figure  6  shows  the  behavior  of  the  OH  stretch  band  as  a  function  of  temperature  based  on  data  from 
Hudgins  et  al.8  As  the  deposition  temperature  increases  from  10K  to  140K,  the  wavelength  of  peak 
absorption  shifts  from  3.02  to  3.06  pm,  a  change  that  has  been  documented  by  many  groups.  The 
absorption  coefficient  increases  by  about  28%.  A  similar  transformation  takes  place  when  ice  is 
deposited  at  25K  and  then  warmed  to  120K  (Raut  et  al.11),  although  there  are  slight  variations. 


Figure  5.  Transmission  spectra  for  a  variety  of  thicknesses  for  ice  deposited  at  80K.  The 
transmissions  greater  than  1 .0  shortward  of  2.5  pm  are  numerical  artifacts. 
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Absorption  coefficient  a  (cm'1) 


Wavelength  (juac) 


Figure  6.  Temperature-dependent  structure  of  the  3-^m  OH  stretch  feature  based  on  tabulations  of 
Hudgins  et  al 8 
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5.  Cubic  Ice 


Cubic  ice  can  form  at  between  137  and  197K,  but  its  optical  constants  are  poorly  known.  Raut  et  al. 
report  a  transmission  spectrum  that  shows  that  the  OH  stretch  feature  splits  into  three  distinct  peaks 
(Figure  7).  They  increase  in  strength  toward  longer  wavelengths  and  lie  at  3.06,  3.14,  and  3.22  pm. 


6.  Evaluating  Signal  in  Arbitrary  Filter  Bands 


The  signal  S  produced  by  a  sensor  can  be  written  in  the  following  simplified  form: 


=  G(a)/ 


P (k)  l(k)  t(X)  dX, 


(V) 


where  the  integral  is  carried  out  over  the  band  of  interest.  G  is  a  gain  function,  P  is  the  source  of 
radiation,  I  is  the  relative  instrumental  response  that  includes  filter  transmission  and  detector 
response,  and  t  is  the  transmission  of  the  ice  film.  In  the  absence  of  ice  (t  =  1 .0),  S  represents  the 
nominal  signal  from  the  sensor.  While  Eq.  (7)  can  be  considerably  more  complicated,  what  we  seek 
here  is  not  an  absolute  response,  but  rather  a  relative  one:  What  is  the  signal  in  the  absence  of  ice 
relative  to  what  it  is  in  the  presence  of  ice? 

By  calculating  the  ratio  71  =  S(with  ice)/S(without  ice),  which  is  equal  to  the  relative  signal ,  we  have 
a  good  measure  of  the  loss:  When  71  =  1 .0,  there  is  no  loss.  When  71  =  0.9,  the  signal  is  90%  of  what 
it  should  be.  When  71=  0.33,  the  signal  is  one  third  of  what  it  should  be.  Figure  8  shows  some 
values  of  71  as  a  function  of  ice  thickness  D  in  three  bands:  3.0-3. 1  pm,  2.8-4. 0  pm  and  11-14  pm. 
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7.  Nuts  and  Bolts  of  Computing  Transmission  Spectra 

Care  must  be  taken  when  using  wavelength  units  that  are  not  expressed  in  cm,  such  as  we  do  here,  i.e. 
pm.  Tabulations  of  a(X)  are  typically  given  in  a(cm_1)  vs  X  (nm).  Regardless  of  the  wavelength 
units  used,  D  must  be  in  inverse  units  of  a;  i.e. ,  if  a  is  in  cm-1,  then  D  must  be  in  cm.  When  the 
optical  constants  are  given,  then  a  can  only  be  correctly  obtained  by  having  D  and  X  in  the  same 
units.  Most  absorption  coefficients  or  optical  constants  are  presented  in  frequency  (wavenumber 
units)  or  do  not  have  equal  intervals  of  wavelength.  For  convenience,  it  is  a  good  idea  to  regrid  the 
spectra  into  equal  wavelength  intervals  before  doing  an  integration  over  a  band. 

When  trying  to  assess  the  relative  signal  in  a  given  band,  the  spectrum  of  the  light  source  must  be 
specified.  This  has  important  consequences  because  a  source  whose  radiance  increases  rapidly  at 
longer  wavelengths  will  give  a  different  result  than  one  that  decreases  at  longer  wavelengths.  This  is 
also  true  of  instrumental  response,  which  includes  any  filters  in  the  system.  Care  must  also  be  taken 
when  applying  transmission  curves  to  radiance.  The  two  spectra  can  be  multiplied  together  without 
worry.  However,  should  there  be  a  desire  to  convert  the  spectra  from  wavelength  to  frequency  units, 
the  shape  of  the  spectrum  and,  in  particular,  its  peak  wavelength  can  change.12 

The  amount  of  absorption  (1  - 1)  is  not  a  linear  function  of  thickness  D  (Eq.  1).  Therefore,  one  cannot 
simply  interpolate  between  two  transmission  curves  in  Figure  4  to  get  the  resulting  thickness.  The 
actually  computation  must  be  done.  This  is  especially  true  when  there  is  significant  absorption,  such 
as  near  the  peak  of  an  absorption  feature.  The  only  exception  occurs  when  the  optical  depth  is  much 
less  than  unity  everywhere  within  the  band.  In  this  case,  (1  - 1)  ~ x,  which  is  proportional  to  D 
because 

t  =  exp(--t)  =  1  -  x  +  t2/2!  -  tV3  ! .  »  1  -  x  when  x  «  1 .  (8) 

With  the  absorption  coefficient  or  optical  constants  in  hand,  it  is  a  simple  matter  to  compute  t  directly. 
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k  (solid)  and  reflectivity  R  (dotted) 


8.  Surface  Reflections 


The  reflectivity,  R,  of  a  layer  of  ice  is  given  by  Fresnel’s  Equations,  which  reduce  to  a  single  equation 
for  normal  incidence 


R=  [(m  -  l)2  +  k2]/[(m  +  1) 2  +  k2] 


(9) 


Figure  9  shows  the  reflectivity,  R,  of  ice  in  a  vacuum  as  a  function  of  wavelength.  The  total 
irradiance  of  a  signal  passing  through  a  layer  of  ice  is  then 


I  =  Io  (1  -  R)  exp(-aD)  (10) 


Wavelength  fjjtm.) 


Figure  9.  Amorphous  Ice  reflectivity,  R,  in  a  vacuum  (at  80K)  and  imaginary  part  of  the 
index  of  refraction  as  a  function  of  wavelength. 
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For  very  thin  layers  of  ice,  reflective  losses  can  be  of  the  same  order  or  greater  than  the  absorptive 
losses.  And,  of  course,  there  are  reflective  losses  at  the  interface  between  the  ice  and  the  underlying 
substrate.  These  can  be  calculated  from  Fresnel’s  Equations  as  well. 

A  somewhat  more  troublesome  effect  occurs  when  the  ice  thickness  is  less  than  the  skin  depth  ,©.  In 
this  case,  the  electric  and  magnetic  fields  of  the  vacuum-ice  and  ice-substrate  couple  with  each  other, 
and  the  resulting  amplitudes  are  more  difficult  to  predict.  And  it  obviously  depends  on  index  of 
refraction,  which  is  a  function  of  wavelength.  A  thin  layer  of  ice  may  act  as  an  anti-reflective 
coating,  or  a  pro-reflective  coating.  If  the  ice  thickness  is  much  less  than  1/4  of  the  wavelength  of  the 
incident  light,  such  effects  are  minimal.  Dealing  with  thin  coatings  like  ice  is  well  studied  in  the  field 
of  multilayer  dielectric  interference  filters,  but  that  is  beyond  the  scope  of  this  report,  primarily 
because  the  underlying  substrates  can  have  a  wide  variety  of  optical  constants.  As  a  general  rule, 
however,  the  transmission  of  ice  on  a  substrate  is  independent  of  the  substrate  material. 
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8.  Conclusions 


We  have  presented  an  overview  of  ice’s  absorption  properties  with  an  emphasis  on  cryogenic  ice  (T  < 
137K)  in  the  infrared.  Absorption  coefficients  and  transmission  spectra  are  shown  along  with 
numerical  examples  of  signal  strength  within  several  bands.  We  also  discuss  some  practical  aspects 
of  computing  transmission  spectra  including  reflective  losses  and  the  influence  of  extremely  thin 
layers  of  ice. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer”  for  national  security  programs,  specializing  in 
advanced  military  space  systems.  The  Corporation's  Laboratory  Operations  supports  the  effective  and  timely 
development  and  operation  of  national  security  systems  through  scientific  research  and  the  application  of 
advanced  technology.  Vital  to  the  success  of  the  Corporation  is  the  technical  staffs  wide-ranging  expertise  and 
its  ability  to  stay  abreast  of  new  technological  developments  and  program  support  issues  associated  with  rapidly 
evolving  space  systems.  Contributing  capabilities  are  provided  by  these  individual  organizations: 

Electronics  and  Photonics  Laboratory:  Microelectronics,  VLSI  reliability,  failure  analysis, 
solid-state  device  physics,  compound  semiconductors,  radiation  effects,  infrared  and  CCD 
detector  devices,  data  storage  and  display  technologies;  lasers  and  electro-optics,  solid-state 
laser  design,  micro-optics,  optical  communications,  and  fiber-optic  sensors;  atomic  frequency 
standards,  applied  laser  spectroscopy,  laser  chemistry,  atmospheric  propagation  and  beam 
control,  LEDAR/LADAR  remote  sensing;  solar  cell  and  array  testing  and  evaluation,  batteiy 
electrochemistry,  battery  testing  and  evaluation. 

Space  Materials  Laboratory:  Evaluation  and  characterizations  of  new  materials  and 
processing  techniques:  metals,  alloys,  ceramics,  polymers,  thin  films,  and  composites; 
development  of  advanced  deposition  processes;  nondestructive  evaluation,  component  failure 
analysis  and  reliability;  structural  mechanics,  fracture  mechanics,  and  stress  corrosion;  analysis 
and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures;  launch  vehicle  fluid 
mechanics,  heat  transfer  and  flight  dynamics;  aerothermodynamics;  chemical  and  electric 
propulsion;  environmental  chemistry;  combustion  processes;  space  environment  effects  on 
materials,  hardening  and  vulnerability  assessment;  contamination,  thermal  and  structural 
control;  lubrication  and  surface  phenomena  Microelectromechanical  systems  (MEMS)  for 
space  applications;  laser  micromachining;  laser-surface  physical  and  chemical  interactions; 
micropropulsion;  micro-  and  nanosatellite  mission  analysis;  intelligent  micro  instruments  for 
monitoring  space  and  launch  system  environments. 

Space  Science  Applications  Laboratory:  Magnetospheric,  auroral  and  cosmic-ray  physics, 
wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and  ionospheric  physics, 
density  and  composition  of  the  upper  atmosphere,  remote  sensing  using  atmospheric  radiation; 
solar  physics,  infrared  astronomy,  infrared  signature  analysis;  infrared  surveillance,  imaging  and 
remote  sensing;  multispectral  and  hyperspectral  sensor  development;  data  analysis  and 
algorithm  development;  applications  of  multispectral  and  hyperspectral  imagery  to  defense,  civil 
space,  commercial,  and  environmental  missions;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  Earth’s  atmosphere,  ionosphere  and  magnetosphere;  effects  of 
electromagnetic  and  particulate  radiations  on  space  systems;  space  instrumentation,  design, 
fabrication  and  test;  environmental  chemistry,  trace  detection;  atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions,  and  radiative  signatures  of 
missile  plumes. 


